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Abstract 

We construct hydrogen atmosphere models for strongly magnetized neutron stars in 
thermodynamic equilibrium, taking into account partial ionization. The presence of 
bound states affects the equation of state, absorption coefficients, and polar izability 
tensor of a strongly magnetized plasma. Therefore the partial ionization influences 
the polarization vectors and opacities of normal electromagnetic waves, and thus the 
spectra of outgoing radiation. Here we review a model suitable for the most typical 
neutron-star atmospheres and focus on the problems that remain to be solved for 
its extension to other atmospheric parameters. 
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1 Introduction 

Thermal emission from neutron stars (NSs) may potentially be used to directly 
measure the NS surface magnetic field, temperature, and composition, achieve 
a more complete understanding of the evolution of the NSs, and constrain 
the properties of matter and physical processes under extreme conditions. It 
was realized long ago that a NS atmosphere model should properly include a 
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strong magnetic field and partial ionization (see, e.g., IPavlov et all Il99.4 for 
an early review). Models of fully ionized NS atmosphere s with strong magnetic 
fields were construc ted by several re search groups (e.g. , IShibanov et"Z[ Il992l : 



ZaneetaL 



2000; Ho and Lai , 20031 and references therein) . The most recent 
papers highlighted the effects that ma y be important for the atmospheres 
of magnetars: the ion cyclotron feature ( Ho and Lai I2OOII : IZane et ah . 2001 ) 
and vacuum polarization effects, including a conversion of the normal modes 
of radiation propagating in the magnetized atmosphere (|Ho and Lai . I2OO3I: 
Lai and HoL >2002. .2003. 1. 



Early considerations of partial ioni zation in the magnetized NS atmospheres 
(lMillei]ll992t iRaiagopal et al.lll997l : also reviewed briefiy bv lZavlin fc Pavlovl 
20021 ) were not very reliable because of oversimplified treatments of atomic 
physics and nonideal plasma effects in strong magnetic fields. At the typical 
NS atmosphere parameters, the effects of thermal motion of the bound species 
are important. In the 1990s, binding energies and radiative transition rates 
with allowance for the motion effects in strong magnetic fields have been cal- 
culated for the H atom ()Potekhinl . ll994tlPotekhin and Pavlovl . Il997l l. Recently 
these ator nic data have been iniplemented in calculations of thermod ynamic 
functions ( Potekhin et al. . 1999t Pote diin and Chabrier . 20031 2004), radia- 
tive opacities (|Potekhin and Chabrieil . I2003L l2004l ). and spectra (|Ho et al 



2003t ) of the partially ionized H atmospheres of t he NSs. Some re sults have 
been presented at the previous Cospar meeting ( Ho et al. . 20043). Now our 
atmosphere model has been complemented by the effects of the bound states on 



the p olarization properties of the strongly magnetized plasma (jPotekhin et al 
2004^ . Below we briefly summarize the results that allow us to calculate real- 



istic X-ray spectra of thermal radiation from hydrogen NS atmospheres with 
magnetic flelds B ~ 10^^ — 10^^ G and effective temperatures T > 10^'^ K, 
and outline the problems that remain to be addressed at other atmospheric 
parameters and compositions. 



2 The atmosphere model 



We u se the equation of state (EOS) for H in strong magnetic flelds (jPotekhin et al 
I999I ) based on the free-energy minimization method, which ensures the ther- 
modynamic consistency and allows one to determine number fractions of chem- 
ical species, required for opacity calculations. The model takes into account 
all available theoretical results on the moving H atoms and nonideal Coulomb 
plasmas in the magnetic flelds. This EOS has been tabulated and employed 
for ca lculation of opacities for astrophysical use ( Potekhin and Chabrier , 20031 
20041) . 



It is well known that under typical conditions (e.g., far from the resonances) 
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Fig. 1. Opacities Kj (Eq. [1]) in the 0-mode (left panels) and X-mode (right pan- 
els) versus photon energy in the hydrogen plasma dX B = 10^^ G, T = 10^'^ K, 
and p = 1 g cm~^, for 9b = 10° (upper panels) and 45° (lower panels). Solid 
lines: a self-consistent calculation for a partially ionized plasma (70% of neutrals); 
dot-dashed lines: the model of full ionization. 

radiation propagates in a magnetized plasma in the form of two normal modes, 
called the extraordinary (X) and the ordinary (O). The opacity in the mode 
j {j =X,0) depends on the photon frequency cu, magnetic field B, density p, 
temperature T, and the angle 9b between the magnetic field and propagation 
direction. It can be written as 



where (a = —1, 0, 1) are the cyclic coordinates of the polarization vectors of 
the normal modes, and the quantities ka (a = —1, 0, 1) do not depend on Ob- 
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Potekhin and Chabriei ( 20031 ) calculated ka{oj) and evaluated K,j(uj,9 



b) using 

the polarization vectors of the nori nal modes in the fully ionized plasma, Such 



a calculation (dubbed "hybrid" in Potekhin et al.l 2004 ) was employed in our 



previous model of par tially ionized hydroKen atmospheres of the NSs with 
strone; mae; netic fields (|Ho et all l20oil2004al) . 



In the new model fjPotekhin et al■l ■ l2004^ ■ we take into account the influence 



of the bound species on the polarization vectors of the normal modes, making 
use of the Kramers-Kronig relation ^ between the imaginary and real parts of 
the plasma polarizability. Thus the calculation of the polarization vectors and 
opacities of the normal modes has become self-consistent. The calculations of 
thermal spectra of the NSs show that such self-consistent treatment is nec- 
essary if the number fraction of the bound states exceeds several percent. In 
Fig. 1 we compare radiative opacities calculated with and without allowance 
for the bound species for one particular set of plasma parameters, typical near 
the radiative surface of a moderately cool neutron star with magnetic field 
B = 10^^ G, for two 9b values. In the case shown in the figure, the neutral 
fraction is 70%, thus the self-consistent treatment of the opacities is important. 



3 Conclusion and unsolved problems 



The constructed atmosphere models allow us to calculate realistic spectra of 
thermal X-ray radiation from H atmospheres of the NSs with 10^^ G < S < 
10^^ G and T > 10^'^ K. Examples of these spectra are presented elsewhere 



( Potekhin et al.l . 12004^ . There remain the following unsolved problems that 



prevent us from making similarly reliable calculations beyond these limits. 

1. Although the H EOS and opacities have been calculated for B u p to 10^^ G 
and implemented in the atmosphere model (|Ho et al.l . lioO^I . l2004at) . the calcu- 
lated spectra a.t B > 10^"^ G depend on the adopted model of mode conversion 
due to the vacuum resonance and on description of propagation of photons 
with frequencies below the plasma frequency. Both these problems have not 
been definitely solved. Their solution is also important for modeling the low- 
frequency (UV and optical) tail of the spectrum. 



2. At lower T or higher B, H atoms re combine in H„ molecules and eventually 
form the condensed phase (jLail 1200 iL and references therein). Corresponding 
quantum-mechanical data are very incomplete. 



^ Previously this relation has been used by Bulik fc Pavlov (|l996l ) for a neutral gas 
of H atoms in a strong magnetic field. 
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3. At 10^ G < -B < 10^^ G, transition rates of the moving H atom have not 
been calculated previously because of their complexity. The first calculation of 
the energy spectrum appropriate in t his range of B has b een published when 
the present paper was in preparation (|Lozovik and Volkov. .2004) . 



4. At present it is not possible to calculate accurate atmospheric spectra at 
B > 10^^ G for chemical elements other than hydrogen, because of the impor- 
tance of the effects of finite nuclear mass in the strong field regime. Apart from 



(Al 


-Huiai and Schme 


chei . 2003a \h and for the He"*" ion 




Bezchastnov et al. 


1998h. 



5. A more rigorous treatment of radiative transfer in the atmosp here requires 
solving the transfer equations of the Stokes parameters (see, e.g., Lai and Hoi 
20031 for the cases of fully ionized atmospheres). However, since the nonorthog- 
onal features of the modes due to neutral species are diminished by the center- 
of-mass motion, the effect is expected to be small. 

Finally, let us note that the atmosphere model presented here, together with a 



mode l of radiation from the condensed magnetic surface (|van Adelsberg et al. 
200J), has been successfully used for fitting th e spectrum of the isolated neu- 



tron star RX J1856.5-3754 (|Ho et al.1 . lixMbh . 
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